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^^ I A review is given of the theoretical ideas concerning the mysterious pentaquark baryons 

Cn , proposed during the first year after its discovery. We focus on the difficulties involved with 

pH ■ the constituent quark models and the discrepancy between the QCD predictions and most 

r-{ ' of the phenomenological models. 
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§1. Introduction 



The year 2003 will be remembered as a renaissance of hadron spectroscopy. 

Indeed, there were many reports of the discovery of new narrow resonances. The list 

includes the celebrated pentaquark baryons/' the new charmonium state X(3872),^' 

CN ' new Ds mesons,^) and even a narrow K bound nucleus.^) Among these the discovery 

^s ■ of the first pentaquark baryon, ©"*", was the most striking event. 

■^ i The exotic baryon 0^ was discovered in a seminal experiment performed by the 

^^ ' LEPS group at the SPring-8 facility in Harima, Japan. ^^ This particle was identified 

r"| ' in the K~^n invariant mass spectrum in the reaction ^n -^ K~ + 0~^ -^ K~ +K~^ +n, 

QhI which was induced by a SPring-8 tagged photon beam of energy up to 2.4 GeV. That 

Q^' search was motivated by a theoretical paper by Diakonov, Petrov and Polyakov 

1^ ■ (DPP),^) in which they predicted a narrow S = +1 baryon. 

The reported mass of ©^ is 1540 it 10 MeV, and the upper limit of its width 
is 25 MeV. As this baryon has strangeness +1, it should contain at least one s, 
S^ . and thus the minimal number of quarks for Q'^ is five. From the charge and the 

H \ strangeness, u'^cPs is a possibility as the content of 0^. It is therefore commonly 

called the 'pentaquark'. 

The discovery of 0~^ inspired many other experiments, and the existence of 0^ 
was soon confirmed by various groups, including ITEP (DIANA), ^^ JLAB (CLAS)'^^ 
and ELS A (SAPHIR).^' Its discovery was followed by the discovery of yet another 
exotic baryon, S , found by the NA49 group at CERN.^) The particle S (s'^cPu) 
is another manifestly exotic baryon, whose decay into ^~7r~ has been observed at 
the mass M = 1.862 GeV with a width T < 18 MeV. 

Since the discovery of 0~^ was reported, numerous papers attempting to under- 
stand its properties and structure have been uploaded to the preprint server. In this 
article, we pick out some of the prominent achievements among them, attempt to 
explain their content, and also point out some problems. We concentrate on 0~^ , 
first summarizing the experimental facts (§2) and then giving an overview of several 
important theoretical ideas (§§3 and 4). In particular, we discuss the predictions of 
the quark models with regard to the pentaquark states. 
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As QCD is the fundamental theory of hadrons, our goal is to understand ©^ 
from the viewpoint of QCD. Several attempts have been made using the QCD sum 
rule and lattice QCD simulations. In §5, we give an overview of the results and their 
interpretation. We give a conclusion in §6. 

§2. Facts 

A series of experiments have confirmed that the mass of O^ is Mq+ ~ 1540 MeV, 
while the upper limit of the width is 9 MeV, although the accumulated data exhibits 
some dispersion of the central value of the mass.^''^ There is no consensus regarding 
the width, although the HERMES group^^^ claims a non-upper-limit width of order 
10 MeV. (Their mass value is somewhat smaller, 1528 MeV.) Because it decays 
into nK~^ by the strong interaction, the conservation rules guarantee that it has 
a strangeness 5 = +1, baryon number i? = +1, and charge Q = +1. Thus, the 
hypercharge \sY = B + S = 2, and the third component of the isospin is 13 = 0. No 
corresponding pK'^ (I^ = +1) state is observed at the same mass. This leads to the 
conclusion that the isospin of 0^ is / = (see Table I). It is obvious that this state 
cannot be composed of only three quarks. It also seems important that no S* = +1 
baryon state has been observed below the NK threshold, and thus this state seems 
to be the ground state. 

As the flavor SU{3) is an approximate symmetry of QCD, it is likely that 0^ 
belongs to an SU{3) multiplet. The simplest SU{3) irreducible representation (IR) 
for 0^ is 10, which is, in fact, unique if we restrict ourselves to the representations 
produced by 5 quarks (Fig. 1). This is the representation obtained from the Skyrmion 
model, and also it is that assumed in most of the quark model descriptions. We 
therefore assume here that 6*^ is a member of 10. 
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Table I. Tlie quantum numbers of &^ . Tlie spin and parity liave not been determined. 
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Fig. 1. (a) The 10 pentaquarks and (b) tlie mixing of the 10 plus 8 representations for the pen- 
taquarks. 
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Where are the other members associated with 10? An important member is the 
penta-quark A^* (/ = 1/2, Y = \), which is the state with the flavor content: 



AT* = [(2/3)ss + {l/?,){uu + dd)] x {uud or ddu). 



(2-1) 



Because this is a nucleon resonance, it is quite plausible that the present PDG ta- 
ble^^' already includes this state. DPP^^ boldly identified a 1/2+ nucleon resonance, 
A^(1710), as one such possibility. It is, however, noted that the S'[/(3) breaking can 
cause this state to be mixed with other representations. Jaffe and Wilczek (JW)^^' 
proposed the ideal mixing of a pentaquark octet with 10, that is, they predicted the 
two nucleon states 



^ud — ("""" + dd) X {uud or ddu) , 
N* = ss X {uud or ddu). 



(2-2) 



JW further identified these two states as the two lowest 1/2+ resonances, iV(1440) 
and A''(1710). It may, however, be possible that they correspond to, for example, the 
l/2~ states A^(1535) and iV(1650). If so, it is possible that 0+ has negative parity. 
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Table II. The local QCD operators for diquarks and pentaquarks. The color is denoted by a, b and 
c. C = i'-f'^y^ is the charge conjugation matrix. O^j^ is the product of A'' and K with isospin 
1 = 0. 

A similar ideal mixing is plausible for U {S = —1) states, while S {S = —2) 
appear in two multiplets, S* {I = 3/2) in TO and S* {I = 1/2) in 8. They do 
not mix with each other because of the isospin symmetry. The effective number of 
the strange quarks, rig, is the key to understanding the 10 structure and its mixing 
with 8. In Table II, we list the flavor components of the pentaquark states in the 
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form of the local QCD operators. We use "diquark" notation, in which two quarks 
are combined into the color 3 (antisymmetric), the flavor 3 (antisymmetric) and the 
spin (antisymmetric) states. It is seen from Eq. ()2-l() that the average number of 
strange quarks is 1 for (9^, and it increases to Us = 4/3 for N* , Ug = 5/3 for Z"*, 
and Ug = 2 for S* in the 10 representation. 

It is easy to see that the maximal breaking of 5'f/(3) leads to the ideal mixing, 
where ns{N*j) = 0, while ns{N*) = 2. Similarly, ns{S^^) = 1 and ns{U*) = 3. 
Thus, after the ideal mixing, the spectrum should appear as A^^^^, (0"'", i^*^), {N* , 
S3, 2) ■^1/2) ^si where the multiplets in the parentheses have the same value of n^ 
and are degenerate to leading order in mg. JW predicted the masses of the ^^2 
and ^1/2 states to be around 1750 MeV, while the original result of DPP predicts a 
much larger mass for ^g/g. 

§3. Theory I: Soliton model 

3.1. Skyrmion 

The search for the S = +1 pentaquark at SPring-8 was motivated by a theoreti- 
cal prediction of DPP.^' That study is based on the Skyrme soliton (Skyrmion) model 
of the baryon,^^' which utilizes a topological soliton solution in the SU{3) x SU{3) 
chiral symmetric nonlinear sigma model. Their model has specific SU{3) breaking 
terms, indicated by the chiral quark model, but the strengths of the individual terms 
of the effective Lagrangian are determined phenomenologically. 

In the Skyrme model, the classical soliton solution is invariant under neither 
SU(3) nor spatial rotations. The canonical form is the SU{2) hedgehog imbedded 
inU £ SU{3): 

U{f) = exp{iX-rF{r)), (3-1) 

where 

A = (Al, A2, A3) = ( Q Q 

The profile function F(r) represents the radial shape of the soliton and satisfies 
the boundary conditions F{r = 0) = vr and F{r -^ 00) = 0, so that the topological 
winding number is equal to 1, which is identified with the baryon number. To recover 
the symmetries and thus to specify the spin and flavor quantum numbers, we need to 
"rotate" the Skyrmion in SU{3) flavor space and project the appropriate quantum 
numbers. This method of quantization is called the 'rigid rotor quantization'. In 
the SU{3) Skyrme model with Nc = 3, the lowest-lying states are the flavor octet 
and decuplet states, which correspond to the observed ground state baryons^^-* and 
are also consistent with the quark model prediction with three quarks in the baryon. 
The next possible representation, 10, is not composed of three quarks and therefore 
is exotic. Its spin and parity are restricted to 1/2"*" in this model. 

DPP assumed that the known nucleon resonance A^*(1710) is the Y = 1 member 
of the 10 multiplet and derived the (Gell-Mann-Okubo) mass formula 

M = (1890 - 180 Y) MeV, (3-2) 
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using the parameters of the Hamiltonian determined so as to reproduce the octet 
and decuplet baryons. With this formula, the mass of the Y = 2 member [B = 1 
and S = +1) is predicted to be M = 1530 MeV. This member was first called Z* 
but it is now called 0"*". They also calculated the width of 0"*" in the Skyrme model, 
assuming the simplest meson baryon coupling operator. They found that its width 
may be as small as 15 MeV.*-* 

It is important to note that the Skyrme model predicts only positive parity 
states, unless non-collective meson fluctuations are introduced around the Skyrmion. 
There is another constraint, namely that the spin is also specified by the SU{3) 
representation. This constraint is imposed by the Wess-Zumino-Witten term; that 
is, it results from the stipulation that the allowed representation contains a member 
with (I = J, y = 1). Thus, the spin should be equal to the isospin of the Y = 1 
member of the specified representation. ^^-^ Hence, J'^ = 1/2^ is the only allowed 
spin for the 10 flavor state. 

There are additional studies based on the Skyrmion model. Their main achieve- 
ment is to take account of the mixings of other representation due to the SU{3) 
breaking and to fit the parameter values accordingly. Yabu and Ando^^^ carried out 
a full study of the representation mixing in the octet and decuplet baryons some 
time ago, and the basic idea is given in their paper. In the case of 10, mixings of 8 
and 27 have been introduced, but the qualitative features do not change as a result 
of the mixings. ^^' Diakonov and Petrov^^' have reexamined their study after the 
discovery of S* , finding that the choices of the parameter values can be changed so 
that a smaller mass of ^1/2 is obtained. 

We would like to point out here that the Skyrme model in general predicts many 
(sharp) resonances, which may not necessarily correspond to real baryons, because 
the model is assumed to represent QCD at large Nc- There is a consensus that there 
is a significant 1/Nc correction, especially for the excited states. Furthermore, several 
authors have pointed out that the rigid rotor quantization for the 10 representation 
is not consistent with the large Nc expansion. ^'^-''^"'^^ 

3.2. K^ Skyrmion bound state approach 

Itzhaki et al.^-*^) studied an alternative quantization of the SU{3) Skyrmion, 
the kaon-Skyrmion bound state approach, and they compared their method with 
the rigid rotor quantization. This alternative was originally proposed by Callan and 
Klebanov^^) in the context of describing the strange baryon octet and decuplet. They 
pointed out that the WZW term is repulsive (attractive) for K^ iK~), and the K^ 
bound state does not exist in their original work. Itzhaki et al.^^' reexamined this 
case and found that if the kaon mass is as large as 1 GeV, then a bound state may 
appear. In that case, the most probable quantum number is 1/2"'". 

Hosaka^^' pointed out that the chiral bag model exhibits similar qualitative fea- 
tures and predicts a positive parity ground state for 0^ . The chiral bag model has 
a solution of hedgehog form and, in a sense, interpolates between the two extreme 



*' Subsequently, it was pointed out that there is a numerical error in their calculation, and the 
correct value is 30 MeV.^^' 
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pictures, the quark model (MIT bag model) and the Skyrmion model. Hosaka con- 
sidered a configuration of 4 ud quarks and 1 s as a function of the strength of the 
pion field at the surface of the bag. The first three light quarks occupy the lowest 
K = Q state. {K is the grandspin defined hj K = I + J .) The s state is almost 
independent of the pion field. The fourth light quark is found to occupy the negative 
parity K^ = 1~ state if the surrounding pion field is sufficiently strong. This hap- 
pens because the energy of the K'^ = 1~ state goes below the energy of the K'^ = 1^ 
state for a stronger pion field. In this case, the total spin and isospin are 

uudd{I = 0,,r = l~) + s{J'^ = 1/2") -^ ,r = 1/2+ and 3/2+. 

This mechanism is qualitatively consistent with the kaon-Skyrmion bound state. 

The kaon bound state, or the chiral bag model, may be more appropriate than 
the rigid-rotor quantized Skyrmion. But these two approaches cannot be extended 
to all the members of 10. Thus we need new ideas to describe the S* pentaquark. 

§4. Theory II: Quark Models 



As seen above, the predictions obtained from the Skyrmion model apply only 
to states with limited quantum numbers. In contrast, QCD with quarks and gluons 
should have a richer spectrum with various spin and flavor quantum numbers. We 
therefore consider the general quark model in this section and find what it can yield 
for the pentaquark states. 

4.1. Symmetry in quark models 

First we consider the symmetry restrictions on the tetra-quark states composed 
of u'^d'^ quarks. We consider only the 1 = and C = 3 states, so that they are a part 
of the 0+ baryon. With four quarks, this state possesses the following symmetry:*-* 



/ = 



C = 3 



The remaining degrees of freedom corresponds to the spin and orbital motion. 
For simplicity, we consider the classification in the nonrelativistic limit. If we assume 
that the ground state is the orbitally symmetric state with L = 0, then the only 
allowed spin is S" = 1 with [31] symmetry. 

However, the conventional quark model yields the conclusion that the spin de- 
pendent interactions, especially the hyperfine interaction (HF), are important. The 
HF favors lower spin states, and therefore the 5 = 1 tetraquark state is less favored. 
On the other hand, the S = state is allowed only with the mixed orbital symme- 
try, [31], having L = 1. Therefore the combined total angular momentum of the 
tetraquark is j = 1. Such a state is favored by the HF interaction, but the cost is 
its greater kinetic energy. 

*' We here consider the case for Q"*", but the same symmetry classification apphes to other 
pentaquark states. For the other states, we replace / = by the flavor representation, F = 6. 
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Thus, the possible low energy configurations for the 0^ baryon are 

(A) u^(f{L = Q,S = l,f = l+) + s(l/2~) : J^ = 1/2" or 3/2~, 

(B) u^cf{L = l,S = 0,f = l') + s{l/2~) : .P = 1/2+ or 3/2+. 



In either case, there should exist a spin partner, J = 3/2 state, associated with 
the "ground" J = 1/2 state. *•* However, the origin of the J = 1/2 — 3/2 splitting 
differs. For case (A), the spin-spin interaction, oc S{u'^(P) • S{s), is responsible for 
splitting them, while for case (B), the (multibody) spin-orbit force, oc L{u^(P) • S{s), 
is the interaction which causes the splitting. In the ordinary hadron spectrum, the 
spin-spin splitting is much larger than the spin-orbit splitting. Furthermore the 
"three-body" LS force in case (B) is expected to be smaller than the two-body LS. 
According to Dudek and Close, ^^-* the splitting in case (B) is of the order of a few 
tens of MeV, while the spin-spin splitting in case (A) is ~ 200 MeV. In any case, the 
observation of the J = 3/2 partner of 0+ would be critical in oder to distinguish 
among the models of the structure of the pentaquark baryons. 

We consider the above two cases further in order to determine which of the two 
is preferable in the constituent quark model approaches. The main question is which 
of the following two competing energies is larger: 

• Orbital excitation energy: AE{AL = 1) = E{L = 1) - E{L = 0), 

• Hyperfine (HF) interaction: AE{AS = 1) = E{S = 1) - E{S = 0). 

If AE{AL = 1) > AE{AS = 1), case (A) with ,r = 1/2" may give the ground 
state, while if AE{AL = 1) < AE{AS = 1), case (B) with the 1/2+ ground state 
may be realized. 

4.2. Constituent quark model 

A simple constituent quark model is employed here to estimate the above two 
quantities using the knowledge we have regarding the ordinary hadron spectrum. The 
constituent quark is known to be a very useful concept in describing the meson and 
baryon spectra and their properties. It has a mass of order 300 MeV, which is believed 
to be generated by the chiral symmetry breaking in QCD. It can be derived by solving 
the Dyson-Schwinger equation for the quark propagator. ^^^ The interaction kernel 
has strong enhancement in the low momentum region, and it induces dynamical chiral 
symmetry breaking and generates a momentum dependent effective mass. Thus the 
quark acquires the self-energy and is renormalized. It is important to note that the 
conserved currents, such as the isospin, /, the hypercharge, Y, and color charge, C, 
are not renormalized. 

Therefore, the constituent quark has the same quantum numbers as the QCD 
quark but its mass is given roughly by Mb/S. We choose, as typical values of the 
quark mass, m^^^ ~ 360 MeV and nig ~ 540 MeV. Then the sums of the quark 



*' Such a state with the same flavor symmetry does not necessarily exist in the Skyrmion picture. 
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1080 MeV for TV and Z\, 

1260 MeVfor 71, Z'and r*, (4-1) 

1980 MeV for u'^<fs. 



Assuming that the residual interactions are weak,*) these values provide estimates 
of the average masses of the indicated baryons. 

4.3. Orbital excitation 

We first estimate the orbital excitation energy, that is, the energy necessary for 
the excitation of a constituent quark from the lowest single particle energy level to 
the first excited level with L = 1. This energy can be estimated from the excitation 
energy of the negative-parity baryon resonance, e.g., A^*(1535) — A^(940), AE ~ 600 
MeV. This is, however, an overestimate, because it contains the difference of the 
hyperfine interactions in the two baryons. 

A more reliable estimate may be obtained by calculating the difference of the 
kinetic energies in the quark-diquark system. For quarks of constituent mass ~ 360 
MeV, the reduced mass becomes 

/x = 360 X - = 270 MeV. (4-2) 

Assuming that the distance between the quark and the diquark is about 0.6 fm, which 
is a typical size of the quark wave function in the hadron, we obtain AE{AL = 1) ~ 
400 MeV. 

It is interesting that the JW diquark model may require even a larger kinetic 
energy excess, because, if we assume that the diquark mass is about 420 MeV, then 
the reduced mass for the diquarks is 210 MeV. 

4.4. Hyperfine interaction 

The importance of the hyperfine interaction in the hadron spectrum was first 
pointed out by DeRujula, Georgi and Glashow^^^ in the context of one-gluon ex- 
change. It is the spin-spin interaction given by the magnetic component of the 
one-gluon exchange (Fig. 2(a)) and is written in nonrelativistic form as 

-±-{X,.X^){B,.a^)5{nj). (4-3) 

This is known as the color magnetic (CM) interaction. Its main contribution in the 
baryon spectrum is to split the octet and decuplet baryons, such as N and A. 
We simplify the form of the CM as 

^CM = ^ ^CM iij (Aj • \j){c7i ■ aj) (4-4) 

i<j 

*' A confinement force is, of course, required to combine quarks into bound states. This effect is 
assumed to be witfiin tlie constituent quark mass; i.e., m, is regarded as tfie single particle energy 
of the quarks in the ground state. 
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and fix the strength ZicM using the N — A mass difference, which leads to 

Acu = 18.75 MeV. (4-5) 

Here ^ denotes a factor due to the S'[/(3) symmetry breaking, defined by 



^v 



1 for {ij) = (uu), (ud) or (dd), 

rriu/ms for (ij) = (us) or (ds), 

^{mu/nisf for (ij) 



(4-6) 



[ss . 



This is introduced so that the mass splitting (~ 80 MeV) of A and U is reproduced 
by the CM interaction if we take mg/mu ~ 5/3. 




(b) 



a-X 



pion, kaon 



CTlA 



J J 




Fig. 2. (a) CM from the one gluon exchange, (b) FD from the exchange of the pseudoscalar mesons 
and (c) III due to the quark-instanton couphng. 

Other possible origins of the HF interaction have been discussed in the literature. 
One of the approaches predicts pseudoscalar meson exchange interactions between 
quarks inside the baryon (Fig. 2(b)). We call this interaction the "flavor dependent" 
(FD) HF interaction. ^^^ This force depends on the flavor and spin of the quarks and 
is given in simplified form as 



Ufb = -AFr>y^^S,ij{X{ ■ Aj)((Ti • aj 



(4-7) 



Kj 



The strength can be determined again using the N — A mass difference: Z\fd ~ 30 
MeV. 

The third possible origin of the HF interaction is the instanton induced inter- 
action (HI) (Fig. 2(c)). It was shown by Shuryak and Rosner^^^ and Takeuchi and 
Oka^^) that the HI gives the correct HF interaction for the ground-state baryon 
spectrum. Takeuchi^'^^ further showed that the spin-orbit interaction in the P-wave 
baryons is also described well by the HI. The HI is also known as the Kobayashi- 
Maskawa-'t Hooft (KMT) interaction, and it results from the coupling of light quarks 
with the instanton. It is given as follows: 



(ijk) 



1 - ;r (<?i • (Jj + (Jj ■ (Jk+ (Jk- (fi) 



5{fij)5{fjk), 
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mj)- (4-8) 



vll\ = V^^^^Af 



i<j 






It contains the flavor antisymmetric {u — d — s) 3-body repulsion and also the flavor 
antisymmetric 2-body attraction, which is induced by contracting the 3-body force 
by contracting a pair consisting of a quark and an antiquark. 

The two-body part of the III can be simplified into the HF form, 



Ein = AmyA{S.i 



i<j 



1 - ^(<f, 



(4-9) 



and its strength is given by Z\iii ~ 125 MeV. 

Although the III is not uniquely identified in the baryon spectrum, its importance 
is known in the meson spectrum and the dibaryon spectrum, 't Hooft^^^ pointed out 
that the III breaks Ua{^) symmetry and thus gives the r] and r]' mass difference. 
Indeed, the above form of the III was used by Hatsuda and Kunihiro^^' in an NJL 
type model to reproduce the pseudoscalar meson spectrum. The importance of the 
three-body III in the H dibaryon prediction was pointed out by Takeuchi and Oka;^^^ 
that is, the three-body repulsion pushes the H mass above the 2yl threshold. 

Which of the contributions of the HF is most important? It is quite likely that 
the actual HF interaction is a combination of i7cM) -^FD and i^m. But, we would 
like to stress that the axial U{1) symmetry breaking caused by the III manifests itself 
in the pseudoscalar mesons, r/ and i]' . Thus at least a part of the HF interaction 
must come from the III. This fact is also favorable, because the CM interaction can 
be reduced so that an abnormally large ag is not necessary to explain the baryon 
spectrum. 

4.5. 0^ mass 

Now, using the above simplified forms of the HF interactions, we estimate the 
baryon masses and apply the same calculation to the pentaquarks. The masses of 
the nucleon and the A resonance are given in the CM model by 

Mn = Smg + {Ucm)n = 360 X 3 - 150 = 930 MeV, (4-10) 

Ma = Snig + {Ucm)a = 360 X 3 + 150 = 1230 MeV. (4-11) 

Similarly, the SU{3) breaking effects are properly taken into account in that model, 
the A mass is given by 

Ma = 3mq + Am+ {Ucm)a = 360 x 3 -F 180 - 150 = 1110 MeV, (4-12) 

where the mass difference between the s quark and the u and d quarks is given by 

Z\m = m^ - mg ~ 540 - 360 = 180 MeV. (4-13) 

We note again that all the spin dependent interactions are assumed to be included 
in the "single particle" energy, i.e., the constituent quark mass. 
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It is found that the CM is not sufficient for the H dibaryon, a system of six 
quarks, represented by u'^tfs'^, with strangeness —2. The CM predicts a deep bound 
state for H, whose mass is estimated as 

Mh = Gnig + 2Am + {Ucm)h = 360 x 6 + 2 x 180 - 450 = 2070 MeV, (4-14) 

while the AA threshold is 2231 MeV. Since the prediction of the bound H by Jaffe^^' 
in 1977, experimental searches were performed for twenty years with a negative 
result. This is the case in which the 3-body III gives approximately a 160 MeV 
repulsion. •^^•* Although the other baryons are reproduced fairy well by any of the 
three models, the fact that the H dibaryon has not been observed gives reason to 
believe that the instanton is necessary. 

Now we evaluate the mass of the 0~^ using the simple model described above. 
For the negative parity 0^ (l/2~), we obtain 

M0+ = 5mg + Am + {Scu)0+ = 360 x 5 + 180 - 250 ~ 1730 MeV. (4-15) 

This estimate may look too crude because it ignores, for instance, the difference 
between the confinement of the 3-quark baryon and the pentaquarks, and also the 
other spin independent forces. It has been shown, however, that a more realistic 
constituent quark model, such as the MIT bag model, gives similar values for the 
pentaquark masses. ^^-^ Thus the CM model requires an extra attraction of about 200 
MeV for 0+ (1/2"). 

The effect of FD on 0^ has been studied by Stancu and Riska^^' and also by 
Jennings and Maltman'^^^ and Bijker et al.^^^ The latter two works also compare FD 
with CM. It is found there that the FD interaction is slightly less repulsive. 

The instanton plays a crucial role also in 0^ . The III (2-body) is strongly at- 
tractive, while the III (3-body) is slightly repulsive. We find the 0"*" mass employing 
the III to be the surprisingly small value 

Mq+ = bnig + Am+ {Um)0+ = 360 x 5 -F 180 - 510 -F 50 ~ 1520 MeV, (4-16) 

where the large attraction (—510 MeV) comes from the 2-body III, and the last term, 
the 50 MeV repulsion, is the contribution of the 3-body III. Although this result is 
very close to the observed 0^ mass, this estimate may not be fully realistic. This is 
because the strong attraction is estimated perturbatively, assuming that the quark 
wave functions of the 3-quark baryon and the pentaquarks are the same. A realistic 
estimate employing the MIT bag model^^-* has been carried out and shows that the 
0+ mass with the full III is about 1600 MeV. 

In contrast to the situation described above, the positive parity baryon has a 
problem that it typically gives a mass that is 100 MeV larger than that of the negative 
parity state. In the CM case, we have 

Me+ = 5mg + Am + (rcM)©+ + AE{AL = 1) 

= 360 X 5 -F 180 - 620 + 450 ~ 1810 MeV. (4-17) 

The attraction due to the HF interaction is indeed strong, but it is cancelled by 
the kinetic energy of the L = 1 orbital motion. The other models yield the same 
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conclusion for 1/2"^ as long as quark configurations based on the independent particle 
model, or the quark shell model, are considered. 

The above exercises show that (1) the constituent quark model for the l/2~ pen- 
taquark tends to yield a larger mass than the observed O"*", though (2) the instanton 
induced interaction is a promising source of strong attraction for the pentaquark 
and (3) the 1/2"*" state may not realize a lower energy than the 1/2" state in the 
independent-particle quark models. 

There are many model calculations that claim a positive parity state of lower 
energy than l/2~. It is, however, fair to say that the majority of such calculations 
contain an arbitrary parameter that is fitted to the observed ©~^ mass or estimates 
the kinetic energy associated with L = 1 unreasonably small. One such example 
is the estimate in Ref.,3'^)''^°) which employ AE{AL = 1) = 210 - 250 MeV, based 
on the splitting of the A^*(1440) from the A^(940). But it is clearly unjustified to 
assume that the A'^*(1440) is a 2 ^ excited state and AE{AL = 1) = huj. 

One possibility is to deviate from the independent particle picture. Recently, a 
dynamical (but semi-classical) five-quark calculation'^^^ was carried out with a model 
Hamiltonian that contains a linear confinement potential plus the CM interaction. 
That calculation employs the antisymmetrized molecular dynamics (AMD) method, 
and the result shows that the J"" = (l/2+,l/2")/ = and J"" = (3/2")/ = 1 
states are almost degenerate as the lowest energy states. The authors claim that the 
diquark-type correlation causes the 1/2"*" state to have a smaller mass. This is the 
subject of the next section. 

4.6. Diquarks 

In order to decrease the masses of the pentaquark states to the observed values in 
the quark model, Jaffe and Wilczek (JW)^^^ proposed the diquark model. A similar 
model was also proposed by Karliner and Lipkin (KL).^'^) 

The diquark is a strongly correlated quark pair. In QCD, both the gluon ex- 
change interaction and the instanton induced interaction favor the spin-singlet and 
color-antisymmetric diquark combination. 

Mo / = C = 3 J^ = 0+. (4-18) 

The instanton model calculation predicts a [ud]Q diquark mass as small as 420 MeV. 
In the JW model, the structure of the 0~^ is considered in terms of two diquarks 
and a strange antiquark system. It is important that the color-singlet bound state of 
these three objects requires an orbital excitation, either to L = 1 or to a state with a 
nodal radial wave function. This can be understood from symmetry considerations: 
The two [ud]o should be antisymmetric in color so that a total color singlet state is 
obtained, while they should be symmetric overall because they are "bosons." Thus 
the relative motion between the diquarks should have L = 1. The mass of the 0^ is 
estimated as 

M ~ 2Md + ms + AE{AL) ~ 840 + 560 + 450 = 1850 MeV. (4-19) 

In the KL model, the combination of a diquark and a triquark uds is considered. 
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and the O^ mass is estimated. It is also assumed that the relative angular momentum 
is L = 1, but the estimated value of AE(AL) is much smaller, ~ 209 MeV. 

How realistic are the diquarks? An attempt was made in lattice QCD to calculate 
the diquark mass. As the diquark is not color singlet, it is necessary to fix the gauge 
to calculate the mass. Wetzorke and Karsch^^) used the Landau gauge and derived 
the spectrum of the diquarks. A recent calculation employing the maximal entropy 
method for the diquark spectral function yields the conclusion that the mass of the 
[ud]o diquark is about 600 MeV. There is also an approach using the Bethe-Salpeter 
equation within rainbow-ladder QCD^^' that predicts a diquark mass ~ 800 MeV. 

We point out that the diquark picture for the pentaquark system requires reex- 
amination of the nucleon and other ground state baryons with the same footing. An 
interesting quantity is the magnetic moments of the ground state baryons. Because 
the diquark enhancement may break the SU{6) symmetry, the relations among the 
magnetic moments of the octet baryons may not hold. For instance, the ratio of the 
magnetic moments of the proton and neutron is given in the diquark limit by 

fJ'p/fJ'n = eu/ed = -2 (4-20) 

if the quarks do not have anomalous moments. 

In summary, the situation is not at all clear in the "constituent" quark model. 
In fact, the quark model has been tested in the case of the nonexotic mesons and 
baryons, but it may be necessary to modify it to accommodate the pentaquarks and 
other exotic states. The suggested strong correlation realized in the diquarks, if it 
is real, may make it necessary to reconsider the basic dynamics of the constituent 
quark model. The width is another quantity that is to be explained in the quark 
model. The "fall-apart" theory is conventionally employed and predicts a width of 
the pentaquarks that is too large, but the validity of such a theory has not yet been 
confirmed. 

§5. QCD prediction 

In the situation that the quark model cannot provide definite predictions for 
the spin, parity and structure of the pentaquark state, we desperately need a QCD 
calculation. Several attempts have been made to apply QCD directly to this problem. 

5.1. QCD sum rules 

The QCD sum rule*^^-* has been applied to the 0~^ spectrum by Zhu,'^^^ Matheus 
et al.,^^' and Sugiyama et al. (SDO).^^^ The first two approaches consider only the 
chiral even terms of the sum rule and do not consider the parity projection. 

SDO made the first attempt to determine the parity and mass of 0^ from 
the QCD sum rule, assuming that its spin is i. The QCD sum rule is derived by 
computing the relevant correlation function in two ways, with the operator prod- 
uct expansion (OPE) in the deeply Euclidean region, and with a phenomenological 
parameterization of the spectral function. This sum rule relates hadron properties 
directly to the QCD vacuum condensates, such as (qq) and ( — G^), as well as to the 
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other fundamental constants, such as nhg. 

In Table II, possible local operators for the pentaquark systems are listed. We 
employ the following interpolating field operator for O^: 

J{x) = eabcSaS!Cs^{x). (5-1) 

Here, a,b,c,--- are color indices and C = i^"^^^. Also, S^{x) = eabcU^{x)C'^^dij{x) 
and Sc{x) = eabcUa{x)Cdi){x) are the scalar (0+) and the pseudoscalar (0~) ud 
diquark operators. They both belong to color 3 and / = 0. The operator J{x) 
produces a baryon with J = 2, / = and strangeness +1. One of the advantages of 
J{x) is that its coupling to the main continuum state, NK, is believed to be small, 
because it cannot be decomposed into a product of N{3q) and K{qq) operators in 
the nonrelativistic limit. 

It is important that the correlation function of J{x) contains both positive and 
negative parity components. Therefore it is necessary for the parity projection to 
determine the parity of the pentaquark state. Such a technique was developed by 
Jido et al.^^^ and it has been applied to the baryon excited states. Here, we follow 
the same procedure for the pentaquark correlation function. 

We consider the retarded Green's function and choose the rest frame, ^=0, 



i7(go) =y d*xe^«-^z(0|e(x'^) J(x)J(0)|0)|,-.=o. (5-2) 

Then the imaginary part of 11 contains two functions, which are the sum and differ- 
ence of the positive and negative parity spectral functions, p , 

-lmn{qo) = A((7o)7° + Biqo) (5-3) 



IT 



or equivalently. 



A{Qo) = ^{p^{qo) + P~{qo)) 
B{qo) = lip+{qo)-p~{qo)) 



p^{qo)=A{qo)±B{qo). (5-4) 



It happens that the part B in the OPE contains only the chiral-odd terms, such as 
(qq) and {qgscr ■ Gq), and thus the chiral symmetry breaking is responsible for the 
parity splitting. 

The sum rule is obtained by comparing the OPE of the correlation function and 
explicit forms of the spectral functions using analytic continuation. The spectral 
function is commonly parametrized as a pole plus the continuum contribution: 

Pphcn(9o) = |A±|^(5(go - "i±) + ^(^0 - \/sth)PcoNT(9o)- (5-5) 

Here, the continuum part is assumed to be identical to the corresponding OPE 
function above the threshold ^/s^. In order to both enhance the pole part and also 
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suppress the higher-dimensional terms of the OPE, we introduce a weight function 
of the "Borel" type. Then the sum rule is obtained as 

rf^oVF(go)Pphen(9o) = / dqQW{qQ)p%y^{qo). (5-6) 

with W{qo) = exp ( —-tt^ ) , where M is the relevant (Borel) mass scale for the baryon. 

The values of the QCD parameters, ttt.^, (ss), {sgscr • Gs) and (^^G^), are chosen so 
that the strange baryon masses are reproduced in the sum rule. The values we choose 
are as follows: m^ = 0.12 GeV, {ss) = 0.8 x (-0.23 GeV)^, mg = (s^.cr • Gs)/{ss) = 
0.8 GeV^ and (^G^) = (0.33 GeV)l 

The residue of the pole, |A-|-|^, in Eq. H5-5() represents the strength with which the 
IF couples to the physical state, and it should be positive only when the pole is real. 
We use this condition to determine the parity of the pentaquark state. In Fig. 3, we 
plot the OPE side (as a function of M) corresponding to |A-i-pexp(— m^/M^). We 
find that the dimension-5 condensate, {sggCr-Gs), gives a large negative contribution 
to |A+p, which ends up with a nearly zero or even slightly negative value. This 
suggests that the pole in the positive-parity spectral function is spurious. In contrast, 
the large {sgs<J • Gs) contribution makes |A_|^ positive, and we thus conclude that 
the obtained negative-parity state is a real state. 

The mass of the negative parity 0"*" is estimated using the sum rule, and we 
find that the M dependence is rather weak. Therefore the reliability of the sum rule 
seems good. It is, however, sensitive to the choice of the threshold, and the result 
is ambiguous to some extent. We have only confirmed that the result is consistent 
with the observed value. 

Summarizing this part, the QCD sum rule shows that (1) the parity splitting is 
caused by chiral symmetry breaking due to rus, (qq) and (qgsCr-Gq) condensates. (2) 
The parity of 0~^ is determined by the positivity of the correlation function and the 
negative parity state is predicted. (3) The mass of the l/2~0"'" baryon is obtained 
as 1.3 - 1.7 GeV. 

An important issue regarding the QCD sum rule is the choice of the interpo- 
lating field operator. If the IF is inappropriate, the sum rule does not contain the 
pentaquark contribution, but the KN continuum state may dominate. A recent 
study ^^' indicates that the continuum effect on the sum rule might change the con- 
clusion concerning the parity assignment. Further studies on the choice of the IF 
and the sum rule approach are necessary. 

5.2. Lattice QCD 

Several lattice QCD calculations suggest that the positive parity 0~^ is very 
massive. Csikor et al.,^'^^ Sasaki^^^ and Lee et al.^^^ calculated the 0'^ correlation 
function in quenched lattice QCD. They employed different local interpolating field 
(IF) operators that do not contain derivatives. They are independent, but from the 
Fierz transform, it is seen that they are not orthogonal. 

Lee et al. have used the IF operator given by the product of N and K, i.e., 
©^^j^ in Table II. They concluded that the correlator is consistent with the NK 
background and exhibit a resonance state in neither 1/2^ nor 1/2+ . This result 
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Fig. 3. Contributions from the terms of each dimension added up successively for the negative- 
parity and positive-parity sum rules with y/sth = 1-8 GeV. 
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can be understood naturally, because their IF has a large overlap with the NK 
continuum state. 

In contrast, in both Refs.^'^) and,^"*^^ it is claimed that there exists a resonance 
state of l/2~ at a mass consistent with the observed O'^ . Sasaki employed the same 
IF as that used in the QCD sum rule in the previous section, and he found that 
the effective mass extracted from the time dependences of the correlation function 
possesses two plateaus. One corresponds to the NK threshold and the other indicates 
a narrow resonance above the NK threshold. It is important to confirm whether 
such a double plateau structure is not an artifact of lattice QCD. 

The Hungarian group^*^-* used another IF and found that 0^ has J'^ = l/2~ and 
that its mass is consistent with the NK threshold. It is not clear whether they can 



Theoretical Overview of the Pentaquark Baryons 17 

claim a resonance above the threshold or not. 

It is fair to say that the majority of the calculations agree that the ground state 
is characterized by J = 1/2, / = 0, and negative parity, while the positive parity 
state has a mass that is greater at least by a few hundred MeV. These calculations 
also show that the results obtained to this time do not depend strongly on the choice 
of the interpolating field. 

§6. Conclusion 

The conclusions of this article are itemized in the following. 

• QCD predicts the existence of a ©'^ (J'^ = 1/2^, / = 0), although it may not 
be easy to distinguish O^ from the KN threshold in the sum rule or in lattice 
QCD. 

• Many models suggest positive parity baryons. The soliton models provide a 
parametrization of the masses of the "pentaquark" states that is consistent 
with the observed masses. The quark models for the positive parity states 
seem to often use a kinetic excitation energy that is too small, AM = M{L = 
1) - M{L = 0) ~ 200 MeV, while a reasonable estimate is about 400 - 450 
MeV. Unless an exotic idea, such as a strong diquark correlation, is introduced, 
it seems unlikely that the l/2~ state will have a higher energy than the 1/2"'" 
state in the quark model. 

• The quark model predicts a mass of 0"*" that is larger than the observed value by 
about 100 MeV - 200 MeV. The discrepancy is larger for 1/2"*". The idea of the 
strong diquark correlation may require confirmation from the QCD viewpoint. 

• The real challenge is to account for the small width. It was even pointed out 
that consistency with the existing KN phase shift analyses necessitates a tiny 
width, as small as 0.1 MeV, in particular for the negative parity O"*".^^-* A 
mechanism yielding strong suppression of the decay is desperately needed. 

New experiments with better statistics and with different production mecha- 
nisms have already appeared and will continue to appear. Most theoretical models 
predict not one but in fact many pentaquark states, and therefore observations of 
other states are extremely important. Among them, the spin J = 3/2 partner and 
the other members of the antidecuplet, as well as the octet, which may be mixed 
with 10, are particularly anticipated. 

Recently, new theoretical papers on the pentaquark appear daily. We need fur- 
ther study of the width and the production mechanisms as well as realistic (not 
toy-model) calculations based on QCD. Lattice QCD with improved statistics and 
other interpolating field operators is also necessary to confirm the nature of pen- 
taquarks. 

It would also be fascinating to determine whether pentaquark spectroscopy is 
related to other newly-observed narrow hadron states, such as D* and X(3872). It 
would be interesting to study whether the heavy quark plays a key role in these new 
resonances. 
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